Entanglement entropy beyond the free case 
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We present a perturbative method to compute the ground state entanglement entropy for inter- 
acting systems. We apply it to a collective model of mutually interacting spins in a magnetic field. 
At the quantum critical point, the entanglement entropy scales logarithmically with the subsystem 
size, the system size, and the anisotropy parameter. We determine the corresponding scaling pref- 
actors and evaluate the leading finite-size correction to the entropy. Our analytical predictions are 
in perfect agreement with numerical results. 
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In recent years, much effort has been devoted to the 
characterization of ground state entanglement in many- 
particle systems. Especially, its relationship with quan- 
tum phase transitions (QPTs) has been investigated, fol- 
lowingthe seminal works in one-dimensional (ID) sys- 
tems [H, 13, Sj ■ Entanglement is also a key concept for 
quantum information theory Several fundamental 
questions have emerged concerning the universality of 
the observed behaviors, as well as their classification. For 
example, the entanglement entropy is known to scale log- 
arithmically with the subsystem size in ID critical spin 
chains 

BIS Si, though its precise form depends on 
the boundary conditions [9|, Ho| . It is natural to won- 
der how that behavior is modified in higher-dimensional 
systems. 

To address this question, we consider the Lipkin- 
Meshkov-Glick model (LMG) [ll[ of mutually interacting 
spins in a magnetic field, i.e. a system with an infinite co- 
ordination number. Although the model was introduced 
in nuclear physics, it has been used to describe many 
other physical systems such as Bose-Einstein condensates 
[13] or small ferromagnetic particles [HI to cite just a few. 
Its entanglement properties have been analyzed from dif- 
ferent perspectives [13. Hi . Iltl . but its entropy has only 



been studied numerically [17 1. 

The aim of this Letter is to investigate analytically 
the entanglement entropy in the LMG model. First, the 
model is studied in the thermodynamic limit for which 
it can be mapped onto a free bosonic system, allowing 
for an exact evaluation of the entropy [3, Qjl, [2(|. In 
a second step, we address the finite-size corrections to 
the entropy at and away from the transition point. This 
leads us to introduce a perturbative method, because for 
a large but finite number of spins, the bosons are weakly 
interacting. At the critical point, the entropy is found 
to grow logarithmically with the subsystem size as in the 
critical ID XY model 21, 22J1 which is the ID counter- 



part of the LMG model. However, the scaling prefactor 
differs from the ID case and also disagrees with previ- 
ous numerical studies [17j|. We also show that, at the 



quantum critical point, the entropy scales logarithmically 
with the system size and the anisotropy parameter, and 
we compute the associated scaling prefactors. Away from 
criticality, the dependence of the entropy on the subsys- 
tem size is found to differ considerably from the scaling in 
finite-dimensional systems. The validity of our approach 
is confirmed by numerical exact diagonalization results. 

The model — The LMG model describes the collec- 
tive behavior of N spins 1/2 with a mutual anisotropic 
(XY) ferromagnetic interaction, subjected to a trans- 
verse magnetic field h. Introducing the total spin op- 
erators S a = c^/2, where a a are the Pauli matrices, 
and the anisotropy parameter 7, the Hamiltonian reads 



H 



1 

"N 



(Si + 1 S 2 y ) ~hS Z 



ID 



Without loss of generality, we can assume ^ 7 < 1 
and ft. ^ 0. As discussed in the literature, this system 
undergoes a second-order QPT at h = 1, between a sym- 
metric (h > 1) and a broken (h < 1) phase, which is well 
described by a mean-field approach. The corresponding 
classical ground state is fully polarized in the field di- 
rection ((<7*}=1) for ft, > 1, and twofold degenerate with 
(<j\)=h for ft < 1. Entanglement originates from quan- 
tum fluctuations around these classical ground states. 

The entanglement entropy — We wish to analyze the 
entanglement entropy of the ground state \ip), defined as 



£ = -Tr A (p A lnp A ) = -Tr s (p B lnpg) 



(2) 



where p A ,B = T^b,a (IV'KV'I)- This definition relies on a 
splitting of the system into two blocks A and B, of sizes L 
and (N — L) respectively. Here, it amounts to introduce 
ieA,B a a/2' To describe quantum fluctuations 
and thus to compute the entanglement entropy, it is con- 
venient to use the Holstein-Primakoff representation of 
these spin operators [23| 



S? = L/2- 
S A = Via* 



a) a. 



yJl-aia/L=(Sff, 



(3) 
(4) 



2 



(N - L)/2-tfb, 



(5) 



S B = s/W^Ltf yjl - tfb/{N -L) = (Sf ) t , (6) 

with = S^ B ± iS^> B . In this way, the LMG 

Hamiltonian is mapped onto a system of two interact- 
ing bosonic modes a and b. The above transformation 
is valid in the symmetric phase, but can also be used in 
the broken phase, provided one first rotates the z-axis to 
bring it along the classical spin direction [l6| . 

The thermodynamic limit — At fixed r = L/N, the 
Hamiltonian can be expanded in 1/N. At order (1/N)° 
and for h > 1, one gets H = NH^ + + 0(1/N) 
with H [ -^ = -ft/2 and 



H (o) 



1 + 7 2/i-7-l 

4~ + 2 

7-lr 



(a f a + 6 f o) 



r(at 2 + a 2 ) + (l-r)(6t 2 + 6 2 ) 



+2 v / -r(l-r)(a t b t + a6) 



(7) 



This effective bosonic Hamiltonian for the spin excita- 
tions is quadratic, and thus exactly solvable. The reduced 
density matrix can be written as pa = e~ K where, at the 



order we consider here, K reads lg, III 



(o) t 



(o) 



(8) 



The key ingredients leading to this form are: i) the 
eigenvalues of pa are non-negative and smaller than one, 
which explains the exponential form; ii) Wick's theorem 
holds for quadratic Hamiltonians, constraining to be 



quadratic. The three coefficients n 
from the three conditions 



(0) 



can be determined 



TrAPA = 1, ((«M) = («M and ((at )) = (a* >, (9) 

where (O) = (ip\n\i/>) and ((0)> = Tr^(e"^Q). To com- 
pute these expectation values, one simply has to diago- 
nalize and K {a \ Then, the Kj^'s are obtained by 
solving the 3x3 nonlinear system of equations ([9]) . From 
these coefficients, we finally obtain the entropy [20| 



S (o) 



i 



1 



with = ct~ 1/2 J [ret + (1 - r)] [r + a(l - r)] , and 



for ft > 1, 

Al — 7 ' 



(10) 

(11) 
(12) 



It is interesting to note that a = limAr^oo 4(5, 2 )/iV also 
plays an important role for the concurrence which is a 
two-spin entanglement measure 15, 16^ . 

In the broken phase h < 1, the ground state is twofold 
degenerate in the thermodynamic limit. The two quan- 
tum states belong to the maximum spin sector S = N/2 




FIG. 1: (Color online) Entanglement entropy as a function 
of the magnetic field, with fixed 7 = 1/4, r = 1/4 and 
N = 32,64,128,256 (from numerics) and 00 (£ (0) ). Arrows 
indicate the behavior of the finite-size correction in various 
regions. For h < 1, we have plotted + In 2 for the reason 
given in the text. The inset is a zoom around h = y/y where 
£ (0) = 0, for N = 64 (black line) and 00 (red line). 



and are eigenstates of the spin-flip operator FJ i a\ . This 
degeneracy is lifted for finite N. Here, we calculated the 
entropy stemming from quantum fluctuations around one 
of the (fully polarized) classical ground states, which do 
not coincide with the quantum ground states just dis- 
cussed. However, they are closely related and it turns out 
that the difference between £(°> and the actual entropy is 
equal to In 2 j24j]. At h = y^y, the entropy vanishes 
since, there, the two degenerate quantum ground states 
can be chosen as separable [l6| . 

The comparison of £ ^ with numerical results obtained 
from exact diagonalization is shown in Fig. [TJ For large 
N, excellent agreement is observed in both phases. 

Critical scaling of the entropy — The main characteris- 
tic of the entropy £^ is that it is finite for ft 7^ 1 whereas 
it diverges at the critical point, in the vicinity of which 
one has 



£(P) 



-iln|ft-l| + iln[r(l-r)] + iln(l-7) + 



l-xln2 + C(|ft- 1| 1/2 ), 



(13) 



with x = 1 for ft > 1 and x 



< 1. 



The 



5/4 for ft 

dependence on ft differs from the one given in Ref. 
where, numerically, the entropy was found to behave as 
— I In. I ft — 1|. The discrepancy comes from the too small 
investigated system sizes 24{. It is also interesting to 
note that the same prefactor —1/4 was reported in the 
Dicke model 25[ using a different approach and it is very 
likely shared by a large class of collective models. 

As already observed for many physical quantities 
[HI, [1(1, the 1/N expansion of the entropy is singular 
at the critical point h = 1. This is reminiscent of a non- 
trivial scaling behavior that we now discuss. We shall 
use the same scaling hypothesis as in Refs. 15, li| which 



3 



assumes that in the vicinity of the critical point, a physi- 
cal observable $ can be written as the sum of a regular 
and a singular contribution, 



$ N (h, 7 ) = <Z> I ° g (h,j) + <Z> s ™ s (h, 7 ). 



(14) 



Here, singular means that the function and/or its deriva- 
tives with respect to h diverge at the critical point, fol- 
lowing a power law. In addition, one has 



~1, 7) 



(h- 1)^(1 



■a*(0. (is) 



where the scaling variable £ = N(h — 1) 3//2 (1 — 7) -1 ^ 2 is 
introduced. The exponents £^ and n$ are characte- 
ristics of the observables $. Since no divergence can oc- 
cur at finite N, one must have £$(C) ~ C -2 ?*/ 3 and 
consequently 

*™*(fc = i, 7 )^iV-("*+ 2 «*/ 3 )(l-7)« + «*/ 3 . (16) 

To perform such an analysis for the entropy, one is led 
to consider $ = e £ which, contrary to £, behaves as a 
power law at lowest order. Combining Eqs. (jTHJ) - dSJ) , 
one readily identifies = — £J = —1/4, and n$ = 0. In 
the large N limit, one thus predicts 

£(h=l)^XN\TiN + X ^lii(l- 1 )+XrHr(l-r)}, (17) 

with xn = X7 = 1/6 an d Xt = 1/2- This means that, at 
criticality and for fixed N, the entropy scales as ^lnL, 
just like in critical ID systems. Although x 7 is in agree- 
ment with previous numerical results [171 ] . it is not the 
case for Xr which was found to be close to 1/3. 

To check our predictions, we have performed a finite- 
size scaling analysis of these prefactors. As can be seen in 
Fig. [21 the exponents have not yet reached their asymp- 
totic value at N — 2000 which is the largest size analyzed 
in Ref . [ItJ ■ A simple extrapolation of these finite-size re- 
sults to the thermodynamic limit confirms the predicted 
values of XNtX-y an d Xr (dotted lines). 

Our approach similarly predicts xn — 1/6 in the Dicke 
model which is consistent with the value of 0.14 ± 0.01 



obtained numerically 251 ]. 

Finite-size corrections — We shall now check the va- 
lidity of the scaling hypothesis (fT5")) . To this purpose, 
we must at least compute the 1/N correction to the en- 
tropy, which requires to develop a perturbation theory to 
go beyond the free (quadratic) boson case. This calcula- 
tion will be performed in the symmetric phase only, since 
this is sufficient to extract the scaling exponents. The 
approach we have developed in this aim constitutes the 
main contribution of the present work. First of all, one 
has to expands = NH^ + /N + 0{1/N 2 ), 

where H^- 1 ' is quartic in a and b. The idea is then to 
expand K = + K^/N + 0{1/N 2 ), with 



r«„t 



a 1 a + 



4 1} («t 2 



a' a + a' a + n 



2 2 



a' a + 



(18) 
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FIG. 2: (Color online) Exponents x~t an< l Xr as a function 
of 1/N obtained from numerical diagonalization of H. For 
clarity, we plotted 2\ T and 6X7 which are expected to be equal 
to 1 in the thermodynamic limit (dotted lines are guides for 
the eyes). Inset: Entropy as a function of In TV at fixed 7 and 
r. The dotted line has a slope \n = 1/6. 

cannot contain any terms of order higher than four 
in a, because is quadratic and H^- 1 ' quartic in a 

and b. In the framework of diagrammatic perturbation 
theory, all terms correspond to certain vertices. Tracing 
out mode b cannot generate, in , terms of sixth order 
in a as every effective 6-legged vertex originates from the 
contraction of at least two bare 4-legged vertices and is 
thus of order 1/N 2 or higher. As in the quadratic case, 
the Ki can be determined from the conditions 

Tr APA =l, ((oM) = (oM, ((at 2 )) = (at 2 ), (19) 
((at 2 a 2 )) = ( a t 2 a 2 ), ( (a t 3 a)) = (a t 3 a); ((a t 4 )) = (a t 4 ); 

which must be satisfied at order 1/N. These expectation 
values can be evaluated perturbatively 24|. An alter- 
native route which we followed here is to compute them 
by i) diagonalizing the quartic operators + H^/N 
and + KW/N using the canonical transformation 

method described in Ref. [2(| which requires to solve a 
system of 48 linear equations for H and of 6 equations 
for K; ii) solving the resulting linear 6x6 system (fT9)) to 
obtain the coefficients This second step can be done 
numerically, but the full exact solutions of this problem 
cannot be given explicitly. However, our main interest 
being the behavior of the entropy near the critical point, 
we have extracted its leading contribution which reads 

3(1-7) 1/2 

<0[(h- l)- 1 ]. (20) 



8(h- I) 3 / 2 

This correction is in complete agreement with the scaling 
hypothesis (|15|) since, at this order and in the vicinity of 
the critical point, one has 



(21) 



[r(l-r 



ll/2 



1/4 



(1-7) 
(h-1) 1 / 4 



1 - 



3(1-7) 1/2 



8N(h- I) 3 / 2 
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FIG. 3: (Color online) Behavior of iV[£ num - £ (0) ] as a func- 
tion of h for fixed r = 1/4 and TV = 32,64, 128,256 (from 
numerics [black lines]) and oo obtained from the pertur- 
bative expansion [red lines]). Inset: Comparison between 
(dots) and its Taylor expansion (|20[) (solid line). 



As a check for the veracity of this perturbative method, 
we compare in Fig. [3] the numerical and analytical 1 /TV 
corrections to the entropy as functions of h. For in- 
creasing 7V, the numerical corrections converge quickly 
towards the analytical one. The inset is a check of the 
Taylor expansion (|20[) in the vicinity of the critical point. 

Discussion — Away from criticality, the entangle- 
ment entropy of typical finite-dimensional systems is, on 
scales greater than the correlation length, proportional to 
the surface area of the considered subsystem [27J (area 
law). At criticality, logarithmic corrections can occur 
[20I l28l [29l . l3fjj . Non-critical collective models however 
behave differently, due to their infinite coordination num- 
ber. In particular, the expansion of (jTUJ) for /( / 1 and 
small t = L/N yields the scaling £ (°) cx — rlnr + 0{t). 
The scaling of the entropy (fTTj) in the critical collective 
model has to be compared to results obtained for other 
critical spin systems. In ID, the entropy at the critical 
point also scales logarithmically with the subsystem size 
L with a prefactor depending on the universality class Q • 
For example, in the XY model in a transverse field, the 
entropy obeys a scaling law similar to (fl~7|) but with differ- 
ing prefactors 21, 22|. That the ID scaling is also found 
in the critical collective model considered here, may seem 
surprising but is simply due to the fact that the ground 
state reduced density matrix is confined to the S = L/2 
sector of the subsystem Hilbert space [17 1 . 

Finally, let us emphasize that the perturbation theory 
we have developed and applied to the LMG model should 
not be associated to this model. It may be extended to 
more complex systems, provided one can determine 
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